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Abstract 



^ \ I review some of the basic information on the Cosmic Neutrino Background momen- 

tum distribution. In particular, I discuss how present data from several cosmological 
observables such as Big Bang Nucleosynthesis, Cosmic Microwave Background and 
(N ■ Large Scale Structure power spectrum constrain possible deviations from a standard 

Fermi-Dirac thermal distribution. 



> 

g 1 The standard picture 

^ \ The large amount of observations accumulated in the last decades, see as a 

\Q \ highly incomplete list of references, provided an unprecedented improvement in our 

O I understanding of the general features of the observable Universe. In fact, this tremen- 

i-C ■ dous experimental effort, supplemented by the generic predictions of inflationary mod- 

els give a nicely consistent picture which is nowadays customary to refer to as the 

O 

^ . 

c/3 . still unsolved major problems, as the nature of Dark Matter or the so fine tuned 

^ ' value of the cosmological constant, it is remarkable that the overall picture of the 

evolution of the Universe can be described in terms of a relatively small number of 
^ \ parameters defining the ACDM model, namely: the (cold+hot) dark matter den- 

c3 [ sity a;dm = ^dm^^, '^^ the baryon density cuf,, the value of the cosmological constant 

or dark energy density ^a, the normalized value of the Hubble parameter today h, 
the primordial perturbation spectrum tilt Ug and the normalization ln[10^°7^rad] 
where T^rad is the curvature perturbation in the radiation era, the optical depth to 
reionization r, the linear theory amplitude of matter fluctuations at 8 Mpc Ug. 

Of course, more exotic scenarios can be considered by adding other (theoretically 
motivated) free parameters, such as extra relativistic degrees of freedom. Quintessence 
equation of state, non spatially fiat Universe etc, which might help in improving 
the agreement of different experimental observations ^U]. In general, the present 
level of precision of data implies that several ideas, both theoretically motivated or 
suggested by experimental results obtained in different frameworks, can be tested 
by checking their implications at the cosmological level. As an example, the overall 
scale of neutrino masses mo can be bound by comparing data with the expected free 



streaming suppression effect of Large Scale Structure (LSS) power spectrum on scales 
smaller than Inr, the horizon when neutrinos become non-relativistic fUE], 



which gives J2''t^u < 0.68 eV ^Hj, where the sum is over the three neutrino species. 
Notice that this bound is stronger than present constraint from earth-based decay 
experiments and of the order of the sensitivity which will be reached in the near future 



It is well known that in the framework of the Hot Big Bang model we expect 
the Universe to be filled by a large amount of neutrinos, with a Fermi-Dirac dis- 
tribution (in the standard scenario) characterized by a temperature of 1.95 "K and 
density of 112 cm~^/flavor. These relic neutrinos decoupled from the electromag- 
netic plasma quite early in time, when weak interaction rates became slower than 
the Hubble rate for (photon) temperatures in the range 2 -h 4 MeV, just before Big 
Bang Nucleosynthesis (BBN) took place. Unfortunately, the fact that this Cosmic 
Neutrino Background (Cz/B) has today a very small kinetic energy, of the order at 
most of 10~^ eV, and that neutrinos only interact via weak interactions, prevents us 
from any possible direct detection of this background on the Earth (see e.g. |14[ [T3] 
for a recent review). 

Nevertheless, there are several indirect ways to constrain the CuB by looking 
at cosmological observables which are influenced either by the fact that neutrinos 
contribute to the Universe expansion rate at all stages, or also via their interactions 
with the electromagnetic plasma and baryons before their decoupling. In this respect, 
one of the the most sensitive probe is represented by the values of the light nuclide 
abundances produced during BBN. Actually, the final yields of Deuterium, ^Li and 
in particular of ^He strongly depend on the number of neutrino species as well as 
on their distribution in phase space at about 1 MeV when the neutron to proton 
density ratio freezes. In fact, since neutrinos were in chemical equilibrium with the 
electromagnetic plasma till this epoch we know by equilibrium thermodynamics that 
they were distributed according to a Fermi-Dirac function, yet BBN can constrain 
exotic features like the value of their chemical potential [TT)| I17j. 

Once decoupled, neutrinos affect all key cosmological observables which are gov- 
erned by later stages of the evolution of the Universe only via their coupling with 
gravity. A well known example is their contribution to the total relativistic energy 
density, which affects the value of the matter-radiation equality point, which in turn 
influences the Cosmic Microwave Background (CMB) anisotropy spectrum, in partic- 
ular around the scale of the first acoustic peak. Similarly, as already mentioned their 
number density and their masses are key parameters in the small scale suppression 
of the power spectrum of LSS. 

The main question which will be addressed in the following Sections is in my 
opinion particularly intriguing: how present (and future) cosmological observations 
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can prove that indeed neutrinos are thermally distributed? I will first consider the 
non thermal features in the CuB distribution arising from the neutrino decoupling 
stage and then discuss the issue on more general grounds. 



2 Neutrino decoupling 

Shortly after neutrino decoupling the temperature of the electromagnetic plasma 
drops below the electron mass, favoring annihilations that heat the photons. As- 
suming that this entropy transfer did not affect the neutrinos because they were 
already completely decoupled, it is easy to calculate the well-known difference be- 
tween the temperatures of relic photons and neutrinos T/T^, = (11/4)^/'^ ~ 1.40. 
However, the processes of neutrino decoupling and annihilations are sufficiently 
close in time so that some relic interactions between and neutrinos exist. These 
relic processes are more efficient for larger neutrino energies, leading to non-thermal 
distortions in the neutrino spectra and a slightly smaller increase of the comoving 
photon temperature. These distortions have been computed by several authors by 
explicitly solving the related Boltzmann kinetic equations jTHj-|2ni, and result to be 
very small at the level of few percent, so that their direct observation is presently 
out of question. However, they should be included in any calculation of observables 
which are influenced by relic neutrinos. For instance, non-thermal distortions lead to 
an enhanced energy density of relic neutrinos parameterized in terms of the so-called 
effective number of neutrinos j2ZI -^gff 



where is the energy density of photons and pu the total radiation energy density. 
This parameter influences the CMB anisotropics by shifting the matter-radiation 
equivalence, which translates into a different power around the first acoustic peak via 
the Integrated Sachs- Wolfe effect. Future CMB experiments, such as PLANCK j2Hl or 
CMBPOL are foreseen to be sensitive to even a tiny change in the Universe radiation 
content, at the level of percent. 

As a second main effect, the distortion on neutrino distribution affects the pre- 
dictions of BBN. On one hand, the increased value of A^^efr shifts the freezing of 
neutron/proton weak processes, and so the eventual yield of ^He. Furthermore, the 
corresponding thermal averaged cross sections for these processes are also directly 
affected by any change in the electron neutrino distribution function. Both effects 
result into a change of the ^He mass fraction Yp at the level of 10~^ CHI-|21I- which 
is quite small, but it has to be taken into account in precise BBN numerical codes 



Finally, distortions in the neutrino distribution also modifies the present neutrino 
number density, so that the contribution of massive neutrinos to the present energy 
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[misniisi]. 



density of the universe is also changed. 

Recently, the neutrino decoupling stage has been studied taking into account the 
effect of flavor oscillations |H21- In this case the neutrino (antineutrino) ensemble is 
described by 3x3 density matrices 



The diagonal elements correspond to the usual occupation numbers of the different 
flavors, while the off-diagonal terms account for neutrino mixing. 

The equations of motion for the density matrices can be cast in the form [H^ 



where H is the Hubble parameter. The first term in the commutator corresponds to 
vacuum oscillation and is proportional to M^, the mass-squared matrix in the flavor 
basis, related to the diagonal one in the mass basis diag(m^, m^, ml) via the neutrino 
mixing matrix. As a reference, in [321 best-fit values from ref. [Mj were considered 



along with the 3a upper bound on 6*13, sin^6'i3 = 0.047. 

The decoupling of neutrinos takes place at temperatures of the order MeV, when 
neutrinos experience both collisions, described by the collisional integral C [p] and 
refractive effects from the medium. The latter correspond in Eq. ()2.3p to the term 
proportional to the diagonal matrix E, the energy densities of charged leptons. 

In Fig. is shown the evolution of the distortion of the neutrino distribution as a 
function of a; = rrieR for a particular neutrino comoving momentum [y = kR = 10), R 
being the scale factor. At large temperatures or small x, neutrinos are in good thermal 
contact with e^. As x grows weak interactions become less effective in a momentum- 
dependent way, leading to distortions in the neutrino spectra which are larger for u^s 
than for the other flavors. Finally, at larger values of x neutrino decoupling is complete 
and the distortions reach their asymptotic values. Fig. |21 shows the asymptotic values 
of the flavor neutrino distribution, for the cases without oscillations and with non- 
zero mixing. The dependence of the non-thermal distortions in momentum is well 
understood and reflects the fact that more energetic neutrinos were interacting with 

for a longer period. Neutrino oscillations reduce the difference between the flavor 
neutrino distortions, and slightly change the final value of Nes and of the comoving 
photon temperature z = TR. Both effects translate into a small change of the '^He 
mass fraction, as also reported in Table 1. 
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Figure 1: Evolution of the distortion of the Ve and Vx — v^^r spectrum for y — 10. In the case with 
^13 7^ one can distinguish the distortions for (upper line) and Vr (lower line). The line labelled 
with T~f corresponds to the distribution of a neutrino in full thermal contact with the electromagnetic 
plasma. From P2| . 



For completeness, I also report the value of the neutrino distribution function as 
obtained by fitting the numerical results for ^13 = 



feq{y) 



1 + 10"^ 
1 + 10"^ 



(1 - 2.2y + 4:.ly^ - 0.047^ 
-4: + 2.1y + 2Ay^ -0.019 y'^ 
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These expressions can be easily incorporated into the numerical tools such as CMBFAST 
P3] or CAMB [36j used to compute CMB and LSS spectra. In fact, as long as neutrinos 
are still relativistic the net effect of the phase space distribution distortion is only 
via the integrated effect provided by N^s, but a more careful analysis of effects when 
neutrinos become non-relativistic should take into account distortions as a function of 
neutrino momenta. As an example, it is easy to calculate from i\2.5\i the present value 
of neutrino energy density. For the simplest case of (almost) degenerate neutrinos, 
i.e. for the neutrino squared mass scale ttiq much larger than the atmospheric squared 
mass difference Am|i one gets 

2^^m^ (2.6) 
93.14 eV ' ^ ' 

where the value in the denominator is slightly smaller than the analogous result in 
the instantaneous decoupling limit (94.12). 
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Figure 2: Frozen distortions of the flavor neutrino spectra as a function of the conioving momen- 
tum. In the case with ^13 7^ one can distinguish the distortions for (upper hne) and Vt (lower 
hue). From [5^ . 

Table 1: Frozen values of z, the neutrino energy density distortion 5pu^ = Spu^/pu^, N^s 
and AYp including flavor neutrino oscillations. From [32] 
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3 Bounds on non-thermal features in the Cz/B 

Apart from the non thermal features discussed so far, which are expected in the 
framework of Standard Model electroweak interactions, neutrino distribution might 
be quite different if neutrinos also interact with exotic form of matter. Therefore, 
obtaining information on the Cz/B distribution is a way to constrain new physics 
beyond our present knowledge of fundamental interactions. Thus, it is an interesting 
issue to understand to what extent present and future cosmological observations can 
bound such exotic interactions by observing the Cz/B momentum distribution. 

On completely general ground, we can specify the neutrino distribution fa{y) by 
the set of moments Q^^^ 

1 / 4 \ (3+")/3 r 

QL"^ = ^(yy) t^^"- jy'^^fMdy , (3.1) 

where the Q^"^ have been normalized to the standard value of the neutrino tempera- 



ture in the instantaneous decoupling limit = (4/ll)^/^T. Notice that all moments 
can be defined if one assumes that neutrino distribution decays at large comoving 
momentum as exp{—y). This is quite expected since at very high y the shape of the 
distribution is ruled by the behavior imprinted by neutrino decoupling as hot relics 
at the MeV scale. 

If we denote by Pmiy), m being the degree of Pmiv), 



k=0 

the set of polynomials orthonormal with respect to the measure y"^ / {exp{y) + 1) 

ey + 



Pm{y) = Y.^t\' , (3.2) 



/ /i 1 Pn{y)Pm{y) = Snm , (3.3) 



then it follows 



where 



dUy) = ^—Y.F^,mPm{y)dy , (3.4) 



F.,m = T.crQ''V , (3.5) 

i.e. a linear combination of moments up to order m with coefficients c^^^ . 

For a Fermi-Dirac distribution all moments can be expressed in terms of the 
number density Q^^-* or, equivalently, as functions of the only independent parameter 
T,^, the first two moments being related to A'efT and today, respectively 

a;. = fi./i2 = 0.058^^T-3Q(°) , (3.6) 

where we have assumed for simplicity that the three neutrinos share the same distri- 
bution. 

In the following I consider only the first two moments and A'^eff as free and 
independent parameters, to be constrained using cosmological data In fact, the 
only way to decide how many Q^^^ should be included in the analysis can be only 
dictated by the sensitivity of the available observational data to the distortion of the 
neutrino distribution. Presently (and surprisingly), it is already very hard to get 
quite strong constraints on the first two moments. In case future observations would 
reach a higher sensitivity on neutrino distribution, it would be desirable to include 
higher order moments, such as the skewness or the kurtosis, related to Q^^ and Q^^^ , 
respectively. 



For the sake of definiteness in a specific example was worked out in details, 
where neutrinos interact with a light scalar field $ with mass M via the interaction 
lagrangian density 

Cint = ■ (3-8) 

v3 i 

This majoron-inspired model have been also considered in |SH]-|1I] for a very small 
light $ field, even lighter than neutrinos. In this case, neutrino annihilations would 
results into a neutrinoless cosmological model [H] . On the other hand, for m$ > 2mo 
the decays of the unstable $ particles might lead to non trivial features in neutrino 
distribution provided 

i) decays take place out of equilibrium, i.e. for temperatures smaller than the 
decaying particle mass M. 

ii) occur after weak interaction freeze-out, otherwise neutrino interactions with the 
electromagnetic plasma would erase any non thermal features. 

Of course, different scenarios can be considered, as for example the case of unstable 
neutrinos Vh decaying into a (pseudo) scalar particle $ and a lighter neutrino z/;, 

Uh^^yi , (3.9) 

but I notice that at this stage the main issue is rather to understand how present (and 
future) data can constrain the non-thermal contribution to the neutrino background, 
rather than to discriminate among different models. Despite of the fact that in [SZj 
the out of equilibrium $ decay case was taken as the reference model, all results are 
quite general and can be applied to different scenarios as well. 



When the $ particles decay, the neutrino distribution gets an additional contribution, 
which in the narrow width limit and in the instantaneous decay approximation at 
temperature To corresponds to a peaked pulse at y^, = M/(2Tc)) so that 



y^f{y)dy = y" 



ey + 1 



dy 



A 



n 



V2 



exp 



TTCT" 



y-y* 

2a2 



(3.10) 



An example of this non-thermal neutrino spectrum is shown in Figure El The param- 
eter A is given by the comoving $ number density at decay A = n(^{Tj:,)R^ 

Correspondingly, the lower moments expressed in terms of the parameters of Eqs. 
(nTTfl and (jSISl) read 
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(3.11) 
(3.12) 
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Figure 3: Differential number density of relic neutrinos as a function of the conioving momentum 
for the non-thermal spectrum in Equation H3.10|l . The parameters are A — 0.018, y* = 10.5 and 
cr = 1, which corresponds to A'cff — 4. From . 

For sufficiently weakly interacting $ particles, so that they are decoupled from the 
thermal bath since at least the BBN epoch, the largest value of A can be bound by 
BBN as a function of the $ mass and number density. Large values of A and M 
in fact, implies a large value of the $ energy density during BBN, thus affecting the 
Hubble expansion rate and the final abundances of ^He and Deuterium. Using present 
data on these two nuclei, see e.g. [301; this translates into A < 0.1, unless the value 
of M is order MeV or larger, in which case A is much more constrained, see Fig. |3J 

The model described so far can also be tested using CMB and LSS power spectrum. 
The results are summarized in third column of Table |21 For comparison I show also 
the corresponding result for the case where neutrino distribution is assumed to be 
the standard Fermi-Dirac, but the value of N^s corresponding to extra relativistic 
degrees of freedom is taken as a free parameter (second column). In particular, the 
best value of the effective for the three models is shown in the ffist line. Notice that 
both models do not improve significantly the with respect to the standard ACDM, 
despite of their larger parameter space. Notice that the bound on the parameter 
q = co'i,(93.2eV/mo) comes essentially from the BBN prior A < 0.1, since the CMB 
and LSS data alone would be compatible with much larger deviations from a thermal 
phase-space distribution (up to A = 1 at 2-0"). Similarly, is poorly constrained, 
and large values for it are still allowed. In the $ decay scenario this implies that 
these decays can take place in highly out of equilibrium conditions. To -C M, the 
only bound being instead on the scalar particle number density at the BBN epoch. 

The fact that A and are so poorly constrained is due the existence of a degen- 
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Figure 4: The 1 a (thin hnes) and 2 a (thick hnes) BBN bounds on the $ number density 
(normahzed to Tj^) versus mass M in MeV. The regions above the contours would be in disagreement 
with the observed primordial abundances of *Hc or D. From 

eracy direction involving mainly oodm, n^o and A^^efr for the two models in which the 
radiation density near the time of decoupling is a free parameter. Values as large as 
i^dm = 0.23, mo = 1.5 eV or N^s = 9 are still allowed at the 2-a level. An improve- 
ment is expected in removing these degeneracies from future experiments, though not 
particularly dramatic. In Table 3 it is shown the result of a standard Fisher matrix 
analysis assuming a fiducial model as reported in the first line of the Table. The 
forecast is based on the foreseen sensitivity of PLANCK combined with the completed 
SDSS redshift survey with effective volume Ves = 1 h"^ Gpc^ and a free bias. While 
the sensitivity on A^efr is now at the level (or better) than percent, both q and mo are 
not very well constrained, since these two parameters are measurable only from the 
free-streaming effect in the matter power spectrum, while all other parameters have 
a clear signature in the CMB anisotropics. 

4 Conclusions 

The role of Cz/B in cosmology is quite ubiquitous, yet it is quite disappointing that any 
direct detection of background neutrinos will be very hard to achieve. Through their 
influence on several cosmological observables several properties of neutrinos can be 
already constrained, such as their mass, possible finite lifetimes, magnetic moments, 
etc. On the other hand, possible exotic features in their momentum distribution 
would also represent quite a unique imprint of new interactions beyond our present 
understanding of fundamental physics. Unfortunately, present data have not already 
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Table 2: Minimum value of the effective (defined as — 21n£, where C is the 
likelihood function) and la confidence limits for the parameters of the three models 
under consideration. 9 is the ratio of the sound horizon to the angular diameter 
distance multiplied by 100, while j3 is the 2dF redshift-space distortion factor. Finally, 
q = a;,^(93.2 eV/mo). The last five lines refer to parameters which can be obtained by 
the set of independent parameters of the first ten lines. From [TFj . 

reached enough sensitivity to put severe constraints in this respect, despite of the huge 
improvements in the last years, but this sensitivity is likely to be highly improved 
in the near future. Waiting for these new exciting times it is worth scrutinizing new 
theoretical perspectives, keeping in mind as a general warning what Pauli said to his 
friend W. Baade about his neutrino hypothesis ^2] (see D. Haidt in these Proceedings) 

"... Today I have done something which no theoretical physicist should ever do: I 
have predicted something which shall never be detected experimentally..." 
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iVeff 
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fiducial values 


0.96 


0.023 


0.14 
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la error 
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Table 3: Expected errors on the parameters of the ACDM + non thermal neutrino 
model. The first line shows the fiducial values, i.e. the parameter values assumed to 
represent the best fit to the data. The second line give the forecast for the associated 
la errors for PLANCK + SDSS. 
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